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Abstract 

The microscopic spectral density for lattice QCD with two flavors and maximally twisted mass 
O : is computed. The results are given for fixed index of the Dirac operator and include the leading 

order a 2 corrections to the chiral Lagrangian due to the discretization errors. The computation is 
carried out within the framework of Wilson chiral perturbation theory. 
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I. INTRODUCTION 

Large scale numerical simulations of twisted mass lattice QCD [j] are currently investi- 
gated in order to access the deep chiral regime of QCD In the twisted mass formulation 
the standard Hermitian Wilson term is replaced by an anti-Hermitian isospin-violating Wil- 
son term, for a review see Under an axial transformation, this modified Wilson term can 
be tranformed into the standard Wilson term while the mass is transformed into a twisted 
mass term. One advantage of this approach is that the fermion determinant of the two flavor 
Dirac operator is bounded from below due to the twisted quark mass. This offers automatic 
control of the problems which the smallest eigenvalues of the Wilson Dirac operator may 
cause for the numerical stability of standard simulations with Wilson fermions. 

The smallest eigenvalues of the Wilson Dirac operator also play a crucial role for the 
spontaneous breakdown of chiral symmetry 4|, |5[ . Here we compute analytically the density 
of these smallest eigenvalues for twisted mass lattice QCD. This microscopic eigenvalue 
density is uniquely determined by the symmetries of the lattice theory and hence can be 
obtained from the low energy effective theory known as Wilson chiral perturbation theory 
6|-|8|. This effective theory describes the finite volume corrections as well as corrections due 
to discretization errors caused by the nonzero lattice spacing a. To order a 2 , the effects 
of the lattice spacing are parametrized in terms of three additional low energy constants. 
While the value of these constants are specific to the exact implementation on the lattice, 
it is essential to know their values in order to extract physical observables such as the chiral 
condensate, S, and pion decay constant F n . The analytical results for the eigenvalue density 
of the Dirac operator at nonzero twisted mass presented here offer a direct way to test Wilson 
chiral perturbation theory against lattice data. Moreover if the test is successful it provides 
a direct way to measure the additional low energy constants as well as the physical ones. 
Such a test has been carried out for the quenched case with a standard (untwisted) mass in 

nn 

[9L I10l|. Finally, we discuss constraints on the low energy constants from QCD inequalities. 

This paper is organized as follows: To settle the notation the next section gives a brief 
introduction to twisted mass QCD. We then turn to the low energy effective theory, Wil- 
son chiral perturbation theory, in section II III The new results for the microscopic Dirac 
eigenvalue density at maximally twisted mass are presented in section IIVI We discuss the 
constraints on the additional low energy constants from the perspective of QCD inequalites 
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the in sections |Vj Finally we draw conclusions in section I VI I 



II. BASICS OF TWISTED MASS QCD 

Here we briefly recall the basics of twisted mass two flavor QCD in the continuum limit 
as well as on the lattice, see [l| for more details. This also introduces the notation used 
throughout this paper. 



A. Twisted mass in the continuum 

In the continuum formulation the twisted mass fermionic action is given by 

S = J d 4 x diDnJn + m + iz t -f 5 r 3 )ip. (1) 
Under the axial transformation 

i>' = exp(2W7 5 r 3 /2)^, ^' = ^ exp(iuj 5 r 3 /2) (2) 
the mass terms get rotated 

m' = m cos(a;) — z t sin (a;) 

z' t = m sin(a;) + z t cos (a;), (3) 

as follows from 

exp(iw7 5 r 3 ) = cos(w) + 275 r 3 sin(w). (4) 

The continuum covariant derivative term is of course invariant under the axial transformation 
since {75, 7 M } = 0. In the continuum we therefore have 

det(£> M 7 M + m + iz t j 5 T 3 ) = detp^ + ml + i^ff^)- (5) 

Note that the twisted source, z' t , vanishes completely if we make the rotation with tan(w) = 
-zt/m. 

If we simply want to evaluate the partition function at some nonzero mass (and zero 
twisted mass) we could start with both m and z t in the determinant as long as we remember 
that this corresponds to the value 



of the quark mass and zero value of the twisted mass. 
Definition: Maximal twist is obtained at m = with u 



= 7i / 2. For maximal twist 



m = Zt 

z' t = 0, (7) 

so that 

Z{m = 0, z t ; a = 0) = Z(m' = z t , z' t = 0; a = 0). (8) 

B. Twisted mass Wilson fermions on the lattice 

With Wilson fermions on the lattice the discretized covariant derivative 

(XT 

D w = 27m(V m + Vy - y V M V; (9) 

is not ant i- Hermit ian and does not anti-commute with 75. However, D\y is 7,5-Hermitian 

ysDwIs = Dl (10) 

and the product with 75, D 5 (m) = 7 5 (Dvf + m ) is therefore Hermitian. These properties 
are unaltered if one adds a clover term to D\y. 

The main motivation to introduce the twisted mass becomes obvious when we write the 
determinant in terms of the eigenvalues, A|(m), of D^(m) 

det(D w + m + izt^z) = det(D 5 (m) + iz t r 3 ) 

= I1( A K™) + iz t){\){m) - izt) = I1(\ 5 M 2 + (11) 

j j 

The square of the twisted mass sets a lower limit on the terms in the product even when the 
eigenvalues of D 5 (m) are smaller in magnitude than m as happens for a ^ 0. The numerical 
problem with small eigevalues of D 5 is therefore regulated by the twisted mass source. 

Since the Wilson term breaks the axial-symmetry the identification (J6j) for the partition 
function is no longer valid on the lattice. However, as the Wilson term is a cutoff artifact 
one is free to choose the m! as the physical quark mass provided that z' = 0. Therefore, 
if we start with m — 0, it is natural to consider the twisted mass z t as the phyiscal quark 
mass, cf. (JTJ). 
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From Eq. (jTTj) it is then clear that the Dirac spectrum relevant for chiral symmetry 
breaking at maximal twist is that of D^{m = 0) 



d 



log Z{m = 0, izt, —izt] a) 



2z t 



dz t ~°~ v " J ™" A 5 (m = 0) 2 + z 2 

Note that, in the twisted-chiral limit we recover the Banks-Casher 

7rp 5 (A 5 (m = 0) = 0; z t \ a) 



p 5 (A 5 (m = 0),z t ;a). 



4j| relation |36] 



E = lim 



(12) 



(13) 



It is therefore of particular interest to know the analytical form of ps(A 5 (m = 0),z t ]a) in 
the microscopic limit. The microscopic eigenvalue density derived below gives exactly this 
form. 



III. WILSON CHIRAL PERTURBATION THEORY WITH TWISTED MASS 

With the twisted source included the static chiral Lagrangian reads 

VC = Tr(m^U + mU ] ) + Tr{z\nU - z t T Z tf) - Tr(a Wf/ + af/W) (14) 



with the sources 



m = mEV, z t = ZpHV and a = aWsV. 



(15) 



Here we have set W 6 = W 7 = K37(. 



In the microscopic limit (aka. the e-regime) for twisted mass Wilson fermions Uj the 
zero momentum modes of the pion fields factorize from the partition function resulting in 
the ^-dependence 



Z% (m, Zt] a) = / det"(U) e 
h(2) 



vc 



(16) 



Here we have written the expression for a sector with fixed index, u, of the Dirac operator. 
The index is defined through 



sign (k 1 7s I k), 



(17) 



where \k) are the eigenstates of Dw- Note that only the real modes of Dw contribute to 



the index 
D 5 (m 



12| . The index may also be obtained from the flow with m of the eigenvalues of 



IV. THE MICROSCOPIC SPECTRUM WITH TWO MAXIMALLY TWISTED 
FLAVORS 



Here we compute the microscopic spectral density of D 5 (m = 0) relevant for two flavors 
at maximal twisted mass. The computation is carried out for fixed index, u, of the Wilson 
Dirac operator. 

In order to derive this density we employ the graded generating functional with index v. 



This is given by 



14Hl6| 

Z^(Z;a) = ydf/Sdet(zf/)^e + i Trs ^[ c/+c/ " 1 ]) +a2Trg ( c/2+c/ ' 2 ), (18) 
where Z = diag(zz t , — iz t , z, z), and the integration is over Gl(3\l)/U(l). The difference with 



17| is that we now have the twisted mass instead of the standard mass. For a discussion of 



the group manifold we refer to 18]. 



The spectral resolvent is obtained from the graded generating functional by differentiation 
with respect to the z source and a subsequent quench of the additional flavors by the limit 

z z 

G^ n (z, z t \ a) = lim -^-Z^izt, -%z u z, z; a). (19) 

1 z->z az 

Finally, the density of eigenvalues, p§(\ 5 , z t \ a), of D 5 follows from 

Ps(A 5 , z t] a) = (j2^t - ^) ) = \MGI\i{z = -A 5 , z t] a)]^ - (20) 



k I N f =2 



Our main task is therefore to evaluate the graded generating function. In |17| it was 
shown that the generating functional 018|) can be rewritten as 



,2 

feJ=l,2,3 



Z» (7-n) - 6 ^ HsHt ??^} r (l/16a 2 )Trg(S 2 +^ 2 ) -tz/8a 2 J f -is^/Sa? 

^3|i^>«J - (167rfl2)2 ]_J sM B ^ Z) e e dete 

x jdU Sdet(iUYe + ^(su+su^)^ (21) 

where the Berezinian is given by 

R ,qs (jgg - jggKjgg -isi)(is 2 -isi) ( 

&3\1W) = 7- —77 —— — (22) 

(t - lSi)(t - is 2 )(t - IS 3 ) 

and 

S = I 




with s = diag(si, s 2 , s 3 ). 

The integral over U results in the a = generating functional which takes the form 
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20] 



.1 1 



1 



x det 



Z^ 1 (x 1 ,X2,x 3 ,x 4 ;a - 0) - 2- 

/ I v {xi) xfJ„ +2 (^l) X?/„ + 3(xi) 



(24) 



X24 + l(x 2 ) ^2^+2 (#2) xll v+3 {x 2 ) 

x 3 I„ +1 (x 3 ) x 2 3 I u+2 (x 3 ) xll u+3 (x 3 ) 

\ (-iyK u (x 4 ) x 4 (-iy +1 K„ +1 (x 4 ) x 2 4 (-iy+ 2 K u+2 (x 4 ) xl(-iy +3 K v+3 (x 4 ) J 



We can thus write 



zUZ;a) 



dsdt 



(>-<rY 2 J fl 3 |iOZ) 

TTfe(-*Sfc 



5 3|l(^) e (l/16a 2 )Trg(S 2 +^ 2 ) e -t5/8a 2 ^ / -is^/Sa 2 



<k.l=l,2,3 



-t 



^3,i({(4) 1/2 },H 2 ) 1/2 ;a = o) 



(25) 



The next step is to simplify the determinant 



det(e" 



JfcJ=l,2,3 



e .o^j.,^ g ■is 1 Z 2 /8a 2 ^-isxZ-sl&a? 
g-is 2 .Zi/8a 2 g-is 2 -Z2/8a 2 g— is 2 Z i /8a 2 
g-«S32i/8a 2 g-is3^2/8a 2 ^—iszZzj8a 2 



(26) 



Since the other terms in the integrand also combine into an anti-symmetric function of the 
Sfc, all terms in the expansions of the determinant as a sum over permutations give the same 
contributions. In the integrand, we can thus make the replacement 



det(e- iSfc 



is^Zj/fih 2 ' 



fc,j=l,2,3 



y Qg-isiZi/8a 2 -is2Z2/8a 2 — iss,Zi/8a 2 



(27) 



The factor e t ( s i z i+ s 2Z2+s 3 z 3 )/8a 2 j g a^gQ^g^ j n ^ the m i X ed term in the exponent of 
e~i6^ 2 ~ Trg (' s ~ 2 ) . The inverse Berezinian of Z becomes 

1 (z ~ zi)(z - z 2 )(z - z) _ (z - iz t )(z + iz t )(z - z) 

B 3 \i(Z) (z - zt)(z - z 2 )(z 2 - zi) (z - iz t )(z + izt){-2iz t )' 

This contributes a total factor of z/2z 4 to the resolvent G (ie. after differentiation with 
respect to z, and the limit z — > z has been taken, so that we necessarily have to differentiate 
the factor (z — z)). 
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Combining the above expressions the resolvent for D^{m = 0) takes the form 

G u 3ll (z,m = 0,zt;a) (29) 

ds\ds 2 ds 3 dt 



n 2 (16a 2 ) 2 z t Zfj f=2 (iz t , -iz t ] a 

(is 2 - is^jisz - tgiKjgg - is 2 ) 
(t - isi)(t - is 2 )(t - is 3 ) 



x 



(t)» 

x^ii((^) 1/2 ,(^) 1/2 ,(^) 1/2 ) H 2 ) 1/2 ;« = o), 

where the partition function for a = is given by Eq. (|24j) and the integration of S3 + iz is 
over the real axis. The microscopic eigenvalue density of -D5 in the theory with two flavors 
at maximally twisted mass follows from (j20|) . We only need to evaluate the two flavor 
maximally twisted mass partition function which appears in the normalization of Gm. 



A. The microscopic two flavor maximally twisted mass partition function 



In order to complete the computation of the microscopic eigenvalue density we need to 
evaluate the normalization which is given by the two flavor maximally twisted mass partition 
function 



Z%(iz t , -iz t \ a) 



dU det{tUye +1 ^ u+u '^ W ^ u2+u - 2 \ 



(30) 



U{2) 



where Z = dia.g(iz t , —iz-t). Extending the results of 17] to the twisted mass case we find 



Z v 2 {iz u -izt] a) 



le 



la- 



z t Tr(lQa 2 ) 



OO POO 



dsids 2 (isi — is 2 ) 



(31) 



00 j —00 



xe 4« 



K{{si-zt) 2 +(s 2 +zt) 2 } (■ \v(-\v 7 v 



is i y(is 2 yz%(s 1 ,s 2 ;a = 0), 



where 



Z%(x 1: x 2 ; a = 0) 



det 



h{x 2 ) x 2 I u+1 (x 2 ) 



(32) 



The final step in the calculation is to factorize the four dimensional integrals in Eq. ( 1291) 
into the product of two-dimensional integrals. Not only may this factorized form have a 
deep connection to an underlying integrable hierarchy 2l|, it is also highly advantageous 
for numerical evaluation of the eigenvalue density. 
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In Appendix A we show that the spectral resolvent (129]) for the microscopic eigenvalue 
density of -D5 with two flavors at maximally twisted mass can be written as 



G^ x {z,zt\a) = G v 1]{1 {z,z]a) (33) 
Z 2 {iz u z\d) z-iz t riV Z%(-iz t ,z\a) z + iz t 

+ Z»(iz t ,-iz t ;a) 2iz t G ^ l{ %Z * Z] ° j Z»(iz t , -iz t ;a) 2iz t G H^^ 

Here 

G^( Zl ,z 2 ;a) = / dsdt— e -(^)/<i*. a > 

lba z ir J_ QO t + Z2 — ts — z\ 



00 

v 



IS + Z\ 

t + z 2 



Zm(V-(™ + zi) 2 i V~(t + Z2) 2 ,a = 0) (34) 



with 



Z" ll (m 1 ,m 2 ; a = 0) = ( — ] (I u (m 1 )m 2 K u+1 (m2) + m 1 I v+1 {m 1 )K v {m 2 )) (35) 

\m 1 J 

and 

Z 2 ^ 1; z 2 ; a) = — — / d Sl ds 2 -{i Sl + Zl - is 2 - z 2 ) e -^^i)/(i6a 2 ) 

x flSl±zA Z u(^_( i3i + ^ ^J_ {lS2 + Z2? . a = 0) (36) 
\IS 2 + Z 2 J 

with Z 2 (xi,X2;a = 0) given in Eq. (132]) . Note that the first term on the right hand side 
of ( 1331) gives rise to the quenched density of -D5 at zero untwisted mass, m. A similar 
factorization of the unquenched density has been observed in the microscopic limit of QCD 
at nonzero chemical potential j^j . In that case this structure has been understood in terms 
of an underlying integrable hierarchy. 

With Eq. (133]) the spectral density has been expressed in terms of products of double 
integrals. This form is far easier to evaluate numerically than the four fold integral given in 
Eq. ([29]). 

This completes the computation of the microscopic eigenvalue density of D^{m = 0) for 
two flavors at maximal twisted mass in sectors with fixed index of the Wilson Dirac operator. 
See Figured] for plots of the density. Note in particular the behavior of the near zero- modes. 
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^ 5 x 5 

FIG. 1: The spectrum of D§{m = 0) for two flavors with maximally twisted mass. Left: the sector 
with zero index of the Dirac operator. As the twisted mass increases the quenched result (dashed 
curve) is approached. Right: the dependence on the index v for fixed z-t = 4 and a = 0.25. As 
Ws is decreased the near zero-modes become exact (^-functions at A 5 = 0. For small lattice spacing 
the width of the peak is proportional to y/Ws- The thin horizontal line in both plots indicates the 
value 1/ir which is the asymptotic limit of the density for large values of | A5 1 . 

V. QCD INEQUALITIES WITH TWISTED QUARK MASS 

In this section we discuss two QCD inequalities. First, a QCD inequality for the micro- 
scopic partition function in a sector with fixed v and second a QCD inequality for the pion 
masses. We will see that both put constraints on the low energy constants of Wilson chiral 
perturbation theory. 

The twisted mass Nf = 2 QCD partition function is positive definite for all v. This im- 
poses a positivity requirement of the partition function corresponding the chiral Lagrangian 
of the Wilson QCD partition function. Because of the identity 

Zlizt = 0; W 6 , W 7 , W s , a) = (-l) u Z»(z t = 0; -W 6 , -W 7 , -W 8 , a), (37) 

and because for large z t , the sign of the partition function is independent of the Wk, we 
necessarily obtain contraints on the Wk- In case W% = W 7 = we find that W 8 > 0. From 
the small a-expansion of the partition function we obtain the condition 

W 8 - W 6 - W 7 > 0, (38) 



in agreement with the convergence requirements of the graded partition function 15j. Ad- 
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ditional constraints can be obtained from mass inequalities for the pion masses which will 
be discussed in the remainder of this section. 

The Dirac operator including the twisted mass, 

D w + m + i^r 3 75, (39) 

has the Hermiticity property 

Til^{D w + m + ^375)757-1 = (D w + m + ^r 3 7 5 ) t . (40) 
Therefore the inverse Dirac operator 

S(x,y) = (x\ \. \y) (41) 

satisfies 

S(x,yY = 7 5 ti£'(2/,x)tl75. (42) 

Instead of T\ we could of course also have used 

cos(0)ti + sin(0)r 2 (43) 

in the Hermiticity relation ( 14*01 which leads to the same consequences. All we need is a 
combination that anticommutes with T3 and is unitary. This relation allows us to derive 
Weingarten type inequalities 23N25| for the pion masses. 

The correlation function of two messon sources ^>Yi\i{x) and i[jTip{y) evaluated for a fixed 
background gauge field satisfies (r is unitary) 



V i)V rj, i)(x)Ti)(x)i){y)T^(y) = -Tr[S(y, x)TS(x, y)F] + Tr[S(x, x)T]Tr[S(y, y)T] (44) 

= Tr[S'(?/,a;)nri755'(y,x) t zri75r] + Ti[S (x , x)T]Tr[S (y , y)T] 
< Tr[S(y,x)S(y,xy}+Tr[S(x,x)T}Tr[S(y,y)r] . 

The bound in the inequality is saturated for T = 275X1 (or with t\ — > r 2 but not with 
T\ —> r 3 ). This inequality has been evaluated for a fixed gauge field background. However, 
since the fermion determinant is positive for all gauge field configurations the inequality 
continues to hold after averaging. If the disconnected diagrams average to zero, we obtain 

VtfjVtfj tfj(x)TtP(x)ij(y)TtP(y)) < (Tr S(0, x)S(0, x) ] ) . (45) 
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For mesonic channels with mass gap mr we have 

(J V ipV ip ^(x)Tvf)(x)ip(0)Tip(0)) oc exp(— mr|x|) as x — > oo . (46) 

The inequality for the correlators thus translates into an inequality for the meson masses. 
From ( 145]) we then conclude that 26] 



< m-p . 



In particular, we have 



From leading order Wilson chiral perturbation theory one obtains 



27|, 



0\2 



±\2 



16a 2 (W 8 + 2W 6 ) 



(47) 



(48) 



(49) 



If the contribution from disconnected diagrams is not important we conclude that 

W 8 + 2W 6 > 0. (50) 

The contribution of the disconnected diagrams can be isolated by the introduction of valence 
quarks. This results in the inequality [29] 

W 8 > 0, (51) 

independent of the value of W 6 and W 7 . Lattice simulations for twisted mass fermions in 
3(| show that 



mi < rat- 



(52) 



This implies that the contribution of the disconnected diagrams is important for the simula- 



tions in 



301 ] . The possible importance of disconnected diagrams has been studied explicitly 



in lattice simulations of the respective correlators in 



that for the simulations in 



30] 



W 8 + 2W 6 < 0. 



311 ] . Using Eq. ( 149!) we thus conclude 

(53) 



Combined with the inequality (I5TT) derived in 29j we obtain the constraint 

W 6 < 0. 



(54) 



In the quenched case lattice simulations show that the charged pions are the lightest 
pseudoscalar Goldstone bosons 32|. This is an agreement with the lore that disconnected 



diagrams are suppressed in the quenched theory 



25] 
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VI. CONCLUSIONS 



We have computed the microscopic spectral density of the massless Hermitian Wilson 
Dirac operator in the presence of two dynamical flavors at nonzero maximally twisted mass. 
The characteristic shape of the eigenvalue density in sectors with fixed index of the Wilson 
Dirac operator derived in this paper offers a direct way to test Wilson chiral perturbation 
theory for twisted mass against lattice QCD. If the spectral density obtained on the lattice 
follows the analytical prediction, the strong dependence of the analytical result on the low 
energy constant W$ offers a direct way to measure the value of Ws- We have reduced the 
analytical form of the twisted mass microscopic spectral density to a factorized form that is 



easily evaluated with standard numerical methods. A similar 
for two standard dynamical flavors was recently presented in 



actorized form of the density 



33) 



The microscopic results for the spectral density of D 5 at m = have been derived for 



Ws > 0. As has been argued in 
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15] only the theory with W$ > correctly describes 
lattice QCD with Wilson fermions. In support of this we have checked that the microscopic 
partition function for two flavors at maximally twisted mass is a positive definite function 
in all sectors with fixed index v of the Wilson Dirac operator. 
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Appendix A. FACTORIZATION OF G 3{1 

In this Appendix we show that the microscopic eigenvalue density for two flavors of 
maximally twisted mass can be factorized into two-dimensional integrals. 
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We start from the resolvent which is given by Eq. (1291) 



G 3 |iO,m = 0,z t ;a) 
1 



ds 1 ds 2 ds a dt ^- iSl) {) S , 3 -i^r^ (55) 



n 2 {lQa 2 ) 2 Z v Nf=2 {iz u -izt] a) J z t (t - is%)(t - is 2 ){t - is 3 ) 

\(s 1 -z t ) 2 +(s2+zt) 2 + (s3+iz) 2 + (t-z)2]/16a 12 (^l^ggT 



xe 



x^3ii((^) 1/2 ,(^) 1/2 ,(^) 1/2 ,H 2 ) 1/2 ;a = 0) ) 



and use the notation 



{slfl\ fc = l,2,3, 



X4 = it. 



(56) 



Our aim is to rewrite this in a factorized form. To this end we explicitly insert the a = 
partition function given in (}24"1) and consider the combination 

(is 2 - isi)(is 3 - isi)(is 3 - is 2 ) 1 



x det 



(t — is\)(t — is 2 )(t — is 3 ) (£3 — x 2 )(x1 — x\)(x\ — x\) 
( I u {xi) xjv+^xx) x\I u+2 {xi) 

I U (X 2 ) X 2 I u+ i(x 2 ) X%I V+2 {X 2 ) 

I v {x 3 ) x 3 I v+1 (x 3 ) xll v+2 (x 3 ) 



^4+3(^1) 

x 3 2 I u+3 (x 2 ) 
xlh +3 (x 3 ) 



(57) 



V (-\yK u (x 4 ) x A {-\y+ 1 K v+l {x A ) xi(-iy +2 K u+2 (x,) x\(-iy +3 K u+3 (x 4 ) J 

Combining the prefactors and using recursion relations for Bessel functions this can be 
rewritten as 



(xi + x 2 )(x 1 + x 3 ){xx + x 4 ){x 2 + x 3 )(x 2 + x 4 )(x 3 + x 4 ) 
( Iv{xi) xxl v+ i(xi) xll v (xx) 

I U (X 2 ) X 2 I U+1 (X 2 ) X 2 2 l v {x 2 ) 

Iv{x 3 ) x 3 I v+x {x 3 ) x 3 I v {x 3 ) 



(58) 



x det 



x\I v+1 {x 1 ) 
xll u+1 (x 2 ) 
x\I v+1 {x 3 ) 



\{-iyK v {x 4 ) x 4 (-iy +l K u+1 ( X4 ) x 2 (-iyK u (x 4 ) x 3 4 (-iy +3 K u+1 (x 4 ) J 

The factorized form is due to the appearance of I u /K u in the odd columns and I u+ i/K u+1 
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in the even colums. Expanding the determinant results in 

1 



(xi + x 2 )(xi + x 3 )(x 1 + x 4 )(x 2 + x 3 )(x 2 + x 4 )(x 3 + X4) 

X [-(-iy +1 I u+1 (x 3 )K u+ i(x4)Iu(x 1 )I u (x 2 )x 3 X4(xl - x\)(x\ - X$ 

+ (-lYK u+1 (x i )I u+1 (x 1 )I u (x 3 )I u (x 2 )x 1 x 4 (x 2 1 - x\)[x\ - x\) 
-{-l) v K v+1 (x i )I v+ i{x 2 )I v {x 3 )I v {x 1 )x 2 Xi{xl - x\){x{ - x\) 
+ {-l) v+1 I v+ i{x 3 )I v+1 {x 2 )K v (x i )I v {x 1 )x 2 x 3 {xl - xl){x\ - x\) 
-[-\y +1 I vJrl {x 3 )I v j rl {x 1 )K v {x A )I v {x2)x 1 x 3 {x\ - xl)(x 2 2 - x\) 
-{-iyi v j rl {xi)I v+1 {x2)K v (x^I v {x 3 )x 1 x 2 {xl - x 2 2 )(xj - x\)] . (59) 



We then decompose the fractions as 

x 1 x 4y i^x 2 x 3 ) 



(xi + x 2 )(x 1 + x 3 )(xi + x 4 )(x 2 + x 3 )(x 2 + x 4 )(x 3 + X4) 

1 1 



(60) 



(Xi + X 2 )(X 3 + Xi) (Xi + £3) (X 2 + Xi) 

and any cyclic permuatations thereof. This results in 
1 



(Xi + X 3 )(X 2 + Xi) 

x [-(-iy +l I u+1 (x 3 )K p+1 (xi)I u (xi)I p (x 2 )x 3 Xi - (-iyK u+ i(xi)I u+1 (xi)I u (x 3 )I v (x 2 )xiXi 
-(-iy +1 I u+ i(x 3 )I u+1 (x 2 )K u (xi)I u (x 1 )x 2 x 3 - {-iyi v+ x{xx)I v+1 {x 2 )K v {xi)I v {x 3 )x 1 X2\ 

1 



+ 



(X 2 + X 3 )(Xi + Xi) 

x [{-V) vJrl I v+1 {x 3 )K v+1 {xi)I v {x l )I v {x 3 )xiXi + {-iyK v+1 (xi)I I/+ i(x2)I v {x 3 )I u {x 1 )x 2 Xi] 

+ (-l) u+1 I l/+1 (x 3 )I„ +1 (x 1 )K u (xi)I u (x 2 )x 1 x 3 + (-l) i/ ^ + i(a;i)/ l/+ i(x 2 )^(x4)/ 1/ (x 3 )xiX 2 ] 
1 

+ 



(xi + x 2 )(a; 3 + x 4 ) 

x [{-iyK v+1 {xi)I v+ i{xi)I u {x 3 )I v {x 2 )x 1 Xi - (-iyK u+ i(xi)I u+ i(x 2 )I v (x 3 )I u (xi)x 2 Xi 
(-iy +1 I v+1 (x 3 )I v+ i(x 2 )K u (xi)I u (xi)x 2 x 3 - (-l) 1/+1 / i , + i(a; 3 )4 + i(xi)^(x4)4(x 2 )xia; 3 ]. 

_ {-^) v {xiK u+1 {xi)I v {x 2 ) + x 2 A^(x 4 )^ + i(x 2 ) x 3 ^ + i(x 3 )4(xi) - a;i/ y+ i(zi) + h{x 3 ) 

Xi + X 2 Xi + x 3 

(-\y (xjK u+ i(xj)I v (xi) + xiK u (xj)I u+1 (xi) x 2 I v+1 (x 2 )I u (x 3 ) - x 3 I v+1 (x 2 )I v (x 2 ) 
Xi + X\ x 2 + x 3 

(-iy(x 4 K u+1 (xi)I u (x 3 ) + x 3 K„(xj)I u+1 (x 3 ) x 1 I v+1 (x 1 )I v (x 2 ) - x 2 I u+ i{x 2 ) + I u (xx) 

X 3 + Xi X-L + X 2 

(61) 
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Using this identity we can express the resolvent in the factorized form 



G^\ t {z, m — 0, izt, —izt\ a) 



(62) 



G\ 



|i(^ a ^ + ~zf( 



Z%(iz tl z)(z - iz t ) 



G\\ x {-iz u z]d) 



Z 2 {-iz t ,z)(z + iz t ) 



Gi\i(izt,z;a), 



\%z u -iz t )2iz t "~"~'~ v Z 2 (iz t , -iz t )2iz t ~ 1,1 

where Z%(zi, z 2 ) and G^(zi, z 2 \ a) are given in (|36l and ( J34l) respectively. This factorization 
can also be derived in general terms 34j, l35 |. 
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Note that we use the convention of 



151 ] for the low energy constants Wq, W-j and W%. In 
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these constants are denoted by —Wq, —W7 1 and —W% respectively. 
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